We present preliminary measurements of the CP asymmetry parameters 
in KgKgKg decays, reconstructing two of the into tt^tt^ and 

one into In a sample of 227 M BB pairs collected by the BABAR 

detector at the PEP-II B-Factory at SLAC, we find the CP parameters to 
be 5 = -0.25+g:^?(stat) ± 0.05(syst) and C = 0.56+[] j^(stat) ± 0.04(syst). 
Combining this result with the previous BABAR measurement, obtained 
from events with three Kg decaying into tt^tt", we get 

S = -Omtoil (stat) ± 0.04(syst) 
C = -0.10 ±0.25 (stat) ±0.05(syst) 
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1 INTRODUCTION 



In the Standard Model (SM) CP violation arises from the Cabibbo-Kobayashi-Maskawa (CKM) 
quark mixing matrix [1] . Decays of B mesons into charmless hadronic final states with three kaons 
are dominated by 6 ^ sss penguin amplitudes, while other SM amplitudes are suppressed by CKM 
factors [2]. Neglecting these CKM-suppressed contributions, the amplitude of time-dependent 
CP violation for these channels is proportional to sin2/3, where /? = aig{—VcdV*f^/VtdV^l) is the 
CP-violating phase difference between mixing and decay amplitudes and Vij are the elements of 
the CKM matrix. The time-dependent CP-asymmetry is obtained by measuring the proper time 
difference At = tcp — ^tag between a fully reconstructed neutral B meson (Bcp) in the final state 
KgKgKg, and a partially reconstructed recoil B meson (-B^ag)- '^^^ -^tag decay provides evidence 
that it decayed either as B° or B^ (flavor tag). The decay rate f+(f-) when the tagging meson is 
a B'^{B^) is given by 

-\At\/TgO 

[1 it S sin {AnidAt) ^ C COS {AnidAt) ] , 

where r^o is the neutral B meson mean lifetime and Amd is the B^-B^ oscillation frequency. The 
parameters C and S describe the magnitude of CP violation in the decay and in the interference 
between decay and mixing, respectively. The time-dependent CP-violating asymmetry is defined 
asAc-P = (f+-f-)/(f+ + f-). 

Since at first approximation 6 — > s decays can be considered as given by a single amplitude, 
no direct CP violation is expected (C ~ 0) and S ~ —rjf sm2f3, where C (5*) is the parameter for 
direct (mixing-induced) CP violation and r// = -|-1(— 1) corresponds to CP-even (-odd) final states. 
In general, a deviation from these expectations might occur without indicating the presence of 
physics beyond the Standard Model, since a second (CKM suppressed) part is present in the decay 
amplitude. The interference between these two terms can in general produce direct CP violation 
and introduce a nontrivial relation between S and — r/j sin2/3, as a function of the relative phase and 
size of the two amplitudes. It has been noted that for B^ — > K^K^Kg (which is a 77/ = state), as 
for the golden mode S° <^Kg, this suppressed amplitude is a penguin contribution as well [3-5], 
so that the ratio of the two terms is expected to be of the order of A^, where A = 0.2258 it 0.0014 [6] 
is the sine of the Cabibbo angle. 

The value of sin2/3 = 0.726 ± 0.037 determined from tree-level b ccs decays is in good 
agreement with the SM expectaion [6,7]. On the other hand, b sqq processes are dominated by 
one-loop transitions, and hence may have contributions from diagrams with new heavy particles. 
Thus, a sizable deviation of S from — ?7ysin2/3 would be a signal of physics beyond the Standard 
Model [8]. 

Belle and BABAR Collaboration have already reported a measurement of time-dependent CP- 
asymmetry in B^ K^K^^K^g [9]. In the case of BABAR the analysis includes only decaying into 
'K~^'K~ . Since this measurement might be limited in precision by the amount of data one expects 
to have at the end of BABMi experiment, the present work has the main purpose of improving 
the precision using the same dataset, but reconstructing one of the in the 7r°7r° decay mode. 
Because of the absence of charged tracks originating from the B^ decay vertex, we use the vertexing 
technique recently developed for B^ — > K^'k^ [10]. 



8 



2 THE BABAR DETECTOR 



The BABAR detector is described elsewhere [11]. The components are a charged-particle track- 
ing system consisting of a five-layer silicon vertex tracker (SVT) and a 40-layer drift chamber 
(DCH) surrounded by a 1.5-T solenoidal magnet with an instrumented flux return (IFR), an elec- 
tromagnetic calorimeter (EMC) comprised of 6580 CsI(Tl) crystals, and a detector of internally 
reflected Cherenkov light (DIRC) providing excellent charged K — w separation up to a momentum 
of 4.5 GeV/c relevant for this analysis. 

3 ANALYSIS METHOD 

The B^ — > KgKgKg candidate {Bqp) is reconstructed combining three candidates, two of which 
are reconstructed in the tt+tt" mode, while the third is reconstructed in the Kg tt^tt^ 

mode. We reconstruct Kg tt+tt" candidates from pairs of oppositely charged tracks. The two- 
track composites must form a vertex with a tt+tt"" invariant mass within 11 MeV/c^ (about 4a) of 
the nominal Kg mass [12]. We form vr*^ — > 77 candidates from pairs of photon candidates in the 
EMC. Each photon is required to be isolated from any charged tracks, to carry a minimum energy of 
50 MeV, and to have the expected lateral shower shape. We reconstruct Kg candidates from 

7r° pairs which form an invariant mass 480 < m^oTrO < 520 MeV/ c^. B^ — KgKgKg candidates are 
constrained to originate from the e^e~ interaction point using a geometric fit, based on a Kalman 
Filter [13]. We make a requirement on the consistency of the of the fit which retains 93% of the 
signal events, and rejects about 49% of other B decays. We extract the Kg — > tt+tt" decay length 
Lj^o and the invariant mass (m^^) from this fit, and require 100 < m^j < 141 MeV/ and Lj^o 
greater than 5 times its uncertainty. 

For each B candidate we compute two kinematic variables, namely the invariant mass itib and 
the missing mass m,j^jgg = vT^e+e^""^^)^) where q^+e- is the four-momentum of the initial e"'"e~ 
system and qb is the four-momentum of the 5° — KgKgKg candidate after a mass constraint 
on the B^ is applied. By construction the linear correlation coefficient between ^'T^niigg and tub 
vanishes. This combination of variables shows smaller correlation (0.86% on reconstructed signal 
Monte Carlo events and 1.64% on the final data sample) and a better background suppression with 
respect to the equivalent kinematic variables AE and mss used in the BABAR analysis of this mode 
with ah K° vr+vr- in the final state [9]. Using simulations of 5° K^K^K^ and 5° J ji^ir^ 
l^l~)K%{-^ 7r°7r°) decays (/ = e,ii) and reconstructing B^ J/i^{-^ l+l-)K'^{-^ vr'^vr*^) events on 
data, we determine the distribution of ^.u^igg and tub for signal events. We find the signal resolution 
for to be about 40 MeV/ c^, the distribution being asymmetric around the maximum, because of 
leakage effects in the EMC. The signal resolution for JTij^iss' about 6McV/c^, is dominated by the 
beam-energy spread. We select candidates with jub within 150 MeV/ of the nominal B^ mass [12] 
and with 5.11 < m^j^jgg < 5.31GeV/c^. The region m^jgg < 5.2GeV/c^ is devoid of signal and 
used for background characterization. Most background originates from continuum e'^e" — > qq 
{q = u, d, s, c) events, which we suppress using both production and decay properties. To exploit 
the jet-like topology of continuum events, we use the angle 9t between the thrust axis of the Bcp 
candidate and the thrust axis formed from the other charged and neutral particles in the event. 
While I cos^tI is highly peaked near 1 for e'^e" — > qq events, it is nearly uniformly distributed for 
BB events. We require | cos 6*^1 < 0.95. Moreover, we calculate the ratio L2/L0 of two angular 
moments defined as Lj = Y2i \Pi\ \ cos0|p, where p* is the momentum of particle i in the e+e~ rest 
frame, 9^ is the angle between pT and the thrust axis of the B candidate and the sum runs over all 
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reconstructed particles except for the i?-candidate daughters. After all selection requirements are 
applied, the average candidate multiplicity in events with at least one candidate is approximately 
1.67, coming from multiple Kg tt^tt'^ combinations. In these cases, we select the candidate with 
the smallest = ~ ''^k'^)'^ /^nii^ where rrii (mj^o) is the measured (nominal Kg) mass and 

arm estimated uncertainty on the mass of the ith Kg candidate. In simulated events, this 

selection criterion gives the right answer about 81% of the time. The remaining misreconstructed 
events, coming from fake Kg — ^ 7r°7r° candidates, do not affect the determination of At and have 
a small impact on the other variables used in the final fit (the largest correlation is ~ 2.5%). 

Events coming from b ccs would reduce any sensitivity to departures from the Standard 
Model, as this process is characterized by a Standard- Model CP asymmetry (5 ~ sin2/3 and C ~ 0). 
We therefore remove b — *• ccs events by rejecting all candidates with a K^Kg mass combination 
within two times the experimental resolution of the XcO mass. The contribution from Xc2 is found 
to be negligible. Combinatorics from other B decays constitute a further source of background. 
We take this into account by adding a component in the likelihood fit (see Sec. 4), where the shape 
of each likelihood variable is determined from a simulation of inclusive B decays. 

For each B° K^gK^K^g candidate we examine the remaining tracks in the event to determine 
the decay vertex position and the flavor of Stag- Using a neural network based on kinematic and 
particle identification information [14] each event is assigned to one of seven mutually exclusive 
tagging categories, designed to combine flavor tags with similar performance and At resolution. 
We parameterize the performance of this algorithm in a data sample (.Bflav) of fully reconstructed 

average effective tagging efficiency obtained from this sample is 
Q = Ec 4(1 - '^W'f = 0.299 ± 0.005, where and w'^ are the efficiencies and mistag probabilities, 
respectively, for events tagged in category c = 1, 2, • • • 7. For the background, the fraction of events 
(e^) and the asymmetry in the rate of B^ versus 5° tags in each tagging category are extracted 
from a fit to the data. 

The proper-time difference is extracted from the separation of the Bcp and i^tag decay vertices. 
The -B^ag vertex is reconstructed inclusively from the remaining charged particles in the event. To 
reconstruct the Bqp vertex from the single K^ trajectory we exploit the knowledge of the average 
interaction point (IP), which is determined on a run- by-run basis from the spatial distribution 
of vertices from two-track events. We compute At and its uncertainty from a geometric fit to 
the T{4:S) B^B^ system that takes this IP constraint into account. We further improve the 
sensitivity to At by imposing a Gaussian constraint on the sum of the two B decay times (tcp+^tag) 
to be equal to 2r^o with an uncertainty \/2 r^o, which effectively constrains the two vertices to be 
near the ^(45*) line of flight [10]. We have verified in a Monte Carlo simulation that this procedure 
provides an unbiased estimate of At. Details on the vertexing algorithm can be found in Ref. [13]. 

The per-event estimate of the uncertainty on At reflects the strong dependence of the At 
resolution on the K^ flight direction and on the number of SVT layers traversed by the K^ decay 
daughters. In about 97% of the events at least one of the two K^ which decay into tt'^tt~ have 
both pion tracks reconstructed from at least 4 SVT hits, leading to sufficient resolution for the 
time-dependent measurement. The average At resolution in these events is about 1.0 ps. For 
events which fail this criterion or for which a{At) > 2.5 ps or At > 20 ps, the At information is not 
used. However, since C can also be extracted from flavor tagging information alone, these events 
still contribute to the measurement of C. 
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4 MAXIMUM LIKELIHOOD FIT 



We extract the results from unbinned maximum-likelihood fits to the kinematic, event shape L2/Lq^ 
Ai, and flavor tag variables. We maximize the logarithm of an extended likelihood function 

C{S,C,Ns,NB.NsBfsjBjBB,(i) = e-(^^^+^«+^B^) X 

n [^sfselVs {Xi ,yi;S,C) + NsfBe'kPB {Si ,yi;d) + Nj,^f^^elgPj,^{xi , yt ; a)] x (2) 

n [^s{l - fs)e'sVs{x^; C) + Nb{1 - /b)^^^!^.; «) + ^ssi^ " ")] ' 

ieii 

where I {II) is the subset of events with (without) At information. The Nx {X being signal, con- 
tinuum background, or BB background) represent the X component yield, and fx the fraction of 
events with At information. The probabilities Vx {T^x) products of PDFs for each X hypothe- 
ses, evaluated for each event i from the values of Xi = {m^, mj^^jgg, L2/L0, tag, tagging category} 
and yi = {At, a At}- The remaining parameters of the fit are denoted by a. For the B background 
events, the efficiencies and the mistag probabilities and w'^, respectively, for the tagging cate- 
gory c, are fixed to the same values of the signal events. The obscrvables are sufficiently uncorrclatcd 
that we can construct the likelihoods as the products of one-dimensional PDFs. The PDFs for sig- 
nal are parameterized from simulations of signal events. For background PDFs we determine the 
functional form from data in the sideband regions of the other observables where backgrounds 
dominate. Wc include these regions in the fitted sample and simultaneously extract the parameters 
of the background PDFs along with the fit results. All the parameters of BB background PDFs 
are determined using simulated samples of inclusive B decays. All the parameters of the likelihood 
that arc not determined simultaneously with S and C in the final fit are varied according to their 
uncertainties in order to estimate systematic errors. 

The average Az resolution is dominated by the tagging vertex in the event. Thus, we can 
characterize the resolution with a much larger sample of reconstructed B DX decays (-Bg^^y 
sample), which we use as signal parameterization. The amplitudes for the Bcp asymmetries and 
for the ^ga,v A^-vor oscillations are reduced by the same factor due to wrong tags. Both distributions 
are convoluted with a common At resolution function (RF). The RF is parameterized as the sum 
of two Gaussians with a width proportional to the reconstructed a At: and a third Gaussian with 
a fixed width of 8 ps, which accounts for a small fraction of outlying events [14]. The first two 
Gaussians have a non-zero mean, proportional to a At, to account for the small bias in At from 
charm decays on the -Stag ^i*^^- Backgrounds are accounted for by adding terms to the likelihood, 
incorporated with different assumptions about their At evolution and resolution function. 

The fit procedure was tested with both a parameterized simulation of a large number of data- 
sized experiments and a full detector simulation. The likelihood of our data fit agrees with the 
likelihoods from fits to the simulated data. 



5 SYSTEMATIC STUDIES 

Wc obtain systematic uncertainties in the CP coefficients S and C due to the parameterization of 
PDFs for the event yield in signal and background by varying the parameters within one standard 
deviation (evaluated from a fit to Monte Carlo simulated events). We evaluate the uncertainties 
associated with the assumed parameterization of the At resolution function for signal and BB- 
background, a possible difference in the efficiency between B^ and B^, and the fixed values for 
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Am^ and r^o by varying the parameters within one standard deviation (extracted from a fit to the 
-^flav sample). The sum of the two contributions gives the total error associated with the PDF 
parameters. We estimate different uncertainties associated with vertexing. The first is obtained 
by taking the largest value of S'(C)fit — 5'(C)true from fits to signal Monte Carlo events. Here the 
'S'(C)fit represents the result of the fit to our signal Monte Carlo sample, while S'(C)true represents 
input values in the Monte Carlo generation. The second uncertainty is from possible SVT layers 
misalignment. We assign a systematic uncertainty on our knowledge of the beam spot position by 
shifting the beam position in the simulation by ±20 iim in the vertical direction. The sensitivity 
due to any calibration problems or time-dependent effects is evaluated by smearing the beam-spot 
position by an additional dr20/xm in the vertical direction. We include an additional contribution 
from the comparison of the description of the RF between BC vertexing and nominal vertexing in 
the case of J/ipK^ events. Wc estimate also the errors due to the effect of doubly CKM- 

suppressed decays on the tag side [15]. We add these contributions in quadrature to obtain the 
total systematic uncertainty. The summary is reported in Table 1. The largest contribution is 
related to the knowledge of the PDF parameters. For reference, we note that this effect produces 
a systematic error of ±1.7 events on the signal yield. 





AS{+) 


A5(-) 


AC{+) 




PDF parameters 


0.046 


0.039 


0.029 


0.027 


vertexing method 


0.012 


0.012 


0.025 


0.025 


SVT alignment 


0.004 


0.004 


0.008 


0.008 


beam-spot 


0.003 


0.003 


0.005 


0.005 


data/MC RF 


0.006 


0.006 


0.001 


0.001 


doubly-CKM-supprcsscd decays 


0.001 


0.001 


0.011 


0.011 


total errors 


0.049 


0.041 


0.041 


0.039 



Table 1: Summary of systematic uncertainties on S and C. 



6 SUMMARY OF RESULTS 

In the final sample composed of 2748 B candidates we measure 41.0l|| signal events, 2700 ± 56 
continuum background events and 7l^g SS-background events. Assuming the world average for 
the branching ratio ((6.2±0.9)10~^) [9] and the reconstruction efficiency as estimated from a sample 
of simulated signal events, we expected 45 ± 7 events, which is in good agreement with the result. 
We find this preliminary result on CP parameters: 

S = -0.25t|j:^?(stat) ±0.05(syst) 
C = 0.56lgi^(stat)±0.04(syst). 

Fig. 1 shows the background-subtracted distributions of rra^Qjgg and niB for these events, obtained 
using the sPlot weighting technique [16]. Events contribute according to a weight constructed from 
the covariance matrix for the yields (Ns and Nb) and the probability Vs and Vb for the event, 
computed without the use of the variable that is being displayed. The curves represent the signal 
PDFs used in the fit. We combined this result with the previous BABAR measurement, obtained 
using B^ KgKgKg reconstructed from all three Kg decaying into tt+tt" [9]. The combination is 
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Figure 1: Distribution of the event variable m-niiss (l^ft) and niB after reconstruction with the 
weighting technique described in the text. 



obtained through a simultaneous maximum likelihood fit, which takes into account the correlations 
from the common At PDF. The total systematic error is calculated by summing in quadrature the 
uncorrelated sources of errors and taking the largest contribution from the two analyses in the case 
of common sources of background. In this way, we obtain this preliminary result: 

S = -Omtoil (stat) ± 0.04(syst) 
C = -0.10 ±0.25 (stat) ±0.05, (syst) 

In Fig. 2 we show the distributions of signal events, obtained using the sPlot weighting tech- 
nique [16], in the case of K^K^Kg with one reconstructed by vr'^vr'^ mode (left) and for 
the combined fit (right). The superimposed curves represent the results of the fit in the two cases. 

Considering the present uncertainty, this result agrees with Standard Model expectations. A 
future update of this analysis, including new data collected by BABAR, will help to understand if the 
present hint of pattern in the deviation of 6 — > s penguins from the Standard Model predictions [17] 
is a statistical effect or a signal of new physics. 
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Figure 2: Distributions of At for weighted events with -B^ag tagged as (upper plots) or 
(middle plots), and the asymmetry (lower plots). Left plots are for the subsample with all K'^ — > 
TT^vr", right plots are for the combined fit. The points are weighted data and the curves are the 
PDF projections. 
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